Super-soft symmetry energy encountering non-Newtonian gravity in neutron stars 
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Considering the non-Newtonian gravity proposed in the grand unification theories, we show that 
the stabiUty and observed global properties of neutron stars can not rule out the super-soft nuclear 
symmetry energies at supra-saturation densities. The degree of possible violation of the Inverse- 
Square-Law of gravity in neutron stars is estimated using an Equation of State (EOS) of neutron-rich 
nuclear matter consistent with the available terrestrial laboratory data. 
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The density dependence of nuclear symmetry energy 
Esym{p) is an important ingredient for understanding 
many interesting phenomena in astrophysics, cosmology 
[J, 0, i, H and nuclear physics 0, g 0, i, H- How- 
ever, theoretical predictions on the Esym[p) especially 
at supra-saturation densities are currently rather diverse 
B i, 0, [m, [H, [il- Unfortunately, there is no known 
first-principle guiding the high-density behavior of the 
Esym(p) ■ Presently, while many theories, see, e.g., refs. 
[l,El,[li,[ll,[ll,[i3, predict that the Esym{p) increases 

many other models, see, 

M 0, M, m, m, m, 

291 . 130l . I3ll . |32| | , predict that the Egym (p) first increases 
and then decreases above certain supra-saturation densi- 
ties. The Esvm(p) may even become negative at high 
densities 0, i, i, i, [S, [U, [H, [10, HJ. This latter 
kind of symmetry energy functions are generally regarded 
as being soft. Some (e.g., the UV14-I-TNI in [^ and 
group II in [23|) of them can describe very well all ob- 
served properties of neutron stars (NSs). However, the 
super-soft ones the original Gogny-Hartree-Fock 

(GHF) prediction 23| and group III in [2J|) that quickly 



continuously at all densities, 
e.g., refs. jMl, [H, [111, [H, 



drops to zero around three times the saturation density 
Po either can not keep the NSs stable or predict maxi- 
mum NS masses significantly below 1.4 Mq depending 
on the EOS used for symmetric nuclear matter. Given 
the above theoretical situation, experimental indications 
on the high density Egymip) are thus utmost important. 
Very interesti ngly , circumstantial evidence for a super- 
soft Esym{p) |33| was found very recently from analyz- 
ing the FOPI/GSI experimental data on the 7r^/7r+ ra- 
tio in relativistic heavy-ion collisions [s^l within a trans- 
port model [1^ using the MDI (Momentum-Dcpendent- 
Intcraction) EOS [231 . While the symmetric part of 
the MDI EOS is consistent with the existing terrestrial 
nuclear laboratory datajO, the total pressure of NS 
matter obtained using the super-soft Esymip) (which 
is actually the original GHF prediction) preferred by 
the FOPI/GSI data can not keep neutron stars stable. 



Among possibly many important ramifications in astro- 
physics and cosmology, this finding posts immediately a 
serious scientific challenge: how can the NSs be stable 
with such kind of super-soft symmetry energies? In fact, 
this question has been raised repeatedly and the answer 
has been negative long before any experimental indica- 
tion was available. In the literature, the super-soft sym- 
metry energies were often regarded by some people as 
either "unpleasant", see, e.g., [1^, or "unphysical" , see, 
e.g., 0, [H, [131 . These assertions, of course, are all based 
on the assumption that gravity is well understood. How- 
ever, it is really remarkable that gravity, despite being 
the first to be discovered, is actually still considered by 
far the most poorly understood force [H, [s^, [lOl . In fact, 
in pursuit of unifying gravity with the three other funda- 
mental forces, conventional understanding about gravity 
has to be modified due to either the geometrical effect 
of the extra space-time dimensions predicted by string 
theories and/or the exchange of the weakly interacting 
bosons newly proposed in the super-symmetric extension 
of the Standard Model, see, e.g., refs. [4l[,[42| for reviews. 
Consequently, the Inverse-Square-Law (ISL) of gravity is 
expected to be violated. In stable neutron stars at /3 
equilibrium which is determined by the weak and elec- 
tromagnetic interactions, the gravity has to be balanced 
by the strong interaction. Neutron stars are thus a nat- 
ural testing ground of grand unification theories. In this 
Letter, we show that the super-soft Egymip) preferred by 
the FOPI/GSI data can readily keep neutron stars stable 
if the non-Newtonian gravity is considered. 

The deviation from the ISL of gravity can be character- 
ized effectively by adding a Yukawa term to the normal 
gravitational potential [H, [43|, i.e.. 



^r/ N Goomim2 
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where a is a dimensionless strength parameter, A is the 
length scale and Goo is the universal gravitational con- 
stant. Alternatively, the Yukawa term can also be con- 
sidered as due to the putative "fifth force" [HI, [H, [H 
coexisting with grav ity or a non-universal gravitational 
"constant" [U, |4| of G(r) ^ Goo[l + ae-''/^(l + 
r/A)]. In the scalar/vector boson exchange picture. 
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a = / {AttG oo'm'i) and A = (in natural units). 
The , p and are the boson-baryon coupling con- 
stant, the boson and baryon mass, respectively. To re- 
duce gravity from the ISL, the exchange of a vector boson 
is necessary. It is worth noting that a neutral spin-1 vec- 
tor J7-boson has been a favorite candidate. It is very 
weakly coupled to baryons [4^, can mediate the interac- 
tions among Dark Matter (DM) candidates [13, [3 and 
has been used to explain the 511 keV 7-ray observation 
from the galactic bulge [H, [13, [H . 

According to Fujii [52l |. the Yukawa term is simply part 
of the matter system in general relativity. Consequently, 
the Einstein equation remains the same and only the EOS 
is modified. Within the mean-field approximation, the 
extra energy density due to the Yukawa term is [13, \^ 

^UB = 2^y" pi£i)^^piS2)dxidx2 = (2) 
where V is the normalization volume, p is the baryon 



number density and r = \xi—X2\. The corresponding ad- 
dition to the pressure is then Pub = (l ~ ^fl^) ■ 
Assuming a constant boson mass independent of the den- 
sity, one obtains Pub = £ub = k^P^- purposes 
of the present study, it is sufficient to consider neutron 
stars as simply consist of neutrons (n), protons (p) and 
electrons (e). Including the Yukawa term the total pres- 
sure inside neutron stars is P = P_^^^ +-Pub- For the inner 
and outer crusts we use for P^ the EOS of Carriere et al. 
[53! and that of Baym et al. fs^l. respectively. They are 
smoothly connected to the EOS in the core [53. For the 
latter wc use Pnpc{p,5) = p^ [dEo/ dp + dEsym/ dp6'^]^ + 
— d)pEsyiji{p). The value of the isospin asymmetry 
(5 at /3 equilibrium is determined by the chemical equilib- 
rium condition fi^ = fJ-n — fJ-p = 4(S-Esym(p) and the charge 
neutrality requirement pe ~ 5(1 ^ ^)P- The Eo{p) and 
Esymip) obtained consistently within the modified GHF 
approximation are j27i, , respectively, 
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where Pf = h{3TT^^y^^ is the Fermi momentum for sym- 
metric nuclear matter at density p. The coefficients 
Au{x) = -95.98 -a;^ and Ai{x) = -120.57 a;^. 
The values of the parameters are a = 4/3, B = 106.35 
MeV, Ci = -11.70 MeV, C„ = -103.40 McV and 
A = = Pf{po) [131 ■ The resulting symmetric EOS 
contribution dE^/dp to the pressure is consistent with 
that extracted from studying kaon production and nu- 
clear collective flow in rclativistic heavy-ion collisions us- 
ing hadronic transport models assuming no hadron to 
Quark-Gluon Plasma phase transition up to about 5po 
[a, [^ . The parameter x in Eq. [3] was introduced to vary 
the density dependence of the Esy,n{p) without chang- 
ing any property of symmetric nuclear matter and the 
value of Egymipo) = 31 MeV [131 . Shown in the inset 
of Fig.l are two typical Egymip) denoted as MDIxl and 
MDIxO obtained by using x = 1 and x = 0, respectively. 
While the MDIxO Esym{p) increases continuously, the 
MDIxl Esymip) becomes negative above 3po- Only the 



MDIxl Esymip) can reproduce the FOPI/GSI pionpro- 
duction data within the transport model analysis [33| . It 
is seen that the corresponding MDIxl pressure decreases 
with increasing density as shown with the lowest curve 
in Fig. 1. However, the Yukawa term makes the pressure 
grow continuously with increasing density with a value 
of g'^/p^ higher than about 10 GeV"^. 

Shown in Fig. [2] is the mass-radius relation of 
static neutron stars obtained from solving the Tolman- 
Oppcnheimcr-Volkoff (TOV) equation using the MDIxl 
Esymip) and various values for the g^/p^. The result ob- 
tained using the MDIxO without including the Yukawa 
term is included reference [II]- The causality [3| 
and rotational constraint [13| are also shown. The Ke- 
plcrian (mass-shedding) frequency is approximately (sTj 

i^k « l.OS^jl^y \wh{) ^ ^ kHz. So far, the fastest 
pulsar observed is the J1748-2446ad spinning at 716 Hz 
[Sq . Taking 716 Hz as the Keplerian frequency, the M- 
R relation is restricted to the left side of the rotational 
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FIG. 1: (Color online)The inset shows two typical examples 
(MDIxO and MDIxl) of the density dependence of the nuclear 
symmetry energy. The MDIxl (MDIxO) EOS with (without) 
the Yukawa contribution using different values of the / 
in units of GeV~^ are shown. 



3.0 
2.5 

~ 2.0 

o 
5 

1.0 
0.5 





1 1 1 1 

716Hz/ 


< 1 1 1 1 1 1 1 1 1 1 1 1 
MDIxl 




















-MDIxO 




"■"<!50 




- A° 




's 








100 








■■^.75 .. 








50 




-~ y \ 
1 1 1 1 1 


\25^ 


ri~rrr-rt-r-i 1,1," 





10 



15 



20 

R{km) 



25 



30 



35 



FIG. 2: (Color online) The mass-radius relation of static neu- 
tron stars with the MDIxl (MDIxO) EOS with (without) the 
Yukawa contribution. The static neutron star sequences con- 
strained by the rotational frequency 716 Hz of the J1748- 
2446ad [Hy] are taken from Haensel et al. [57| . 

limit. The latter restricts the value of / to less than 
150 GeV"^. It is seen that to produce a neutron star 
with a maximum mass above 1.4 Mq, the j [i^ has to 
be higher than about 50 GeV^^. More specifically, with 
the MDIxl Esymip) and the g^/fi'^ = 50-150 GcV-^, or 
equivalcntly |a|A^ = (2.6 — 7.8) x lO^m^, neutron stars 
can have a maximum mass between 1.4 and 2.5 Mq and 
a corresponding radius between 12 and 18 km. 

For canonical neutron stars of 1.4 M0, the radius is 
quite sensitive to the value used. Thus, besides the 

accurate measurement of neutron star radii, additional 
measurements related to the mass distribution, such as 
the moment of inertia, will be very useful in setting as- 
trophysical constraints on the Esvrn.(p) and g^ j [i^ ■ Ac- 
cording to Lattimcr and Schutz (59[ , at the slow rotation 
limit the moment of inertia can be well approximated as 

4n 



/ « (0.237±0.008)Mi?2 
Shown in Fig. 3 



1 

is the 



/ as a function of M. 



For 



M = 1.4 Mq, the MDIxO without the Yukawa contri- 
bution gives an / no more than 1.8 x 10^^ kg ■ m^ [60| . 
However, significantly larger / values are obtained with 
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FIG. 3: (Color online) The momenta of inertia of neutron 
stars with the MDIxl (MDIxO) Esym[p) with (without) the 
Yukawa contribution. The numbers above the lines are the 
g'^ / values in units of GeV~^. 

the MDIxl Esym{p) and the Yukawa contribution. The 
discovery of the double-pulsar system PSR J0737-3039 
AEzB provides a great opportunity to determine accu- 
rately the moment of inertia I a of the star A [b^, [63j . 
Our results shown here add to the importance of mea- 
suring the moment of inertia accurately. 

To constrain the values of a and A has been a long- 
standing goal of many terrestrial experiments and as- 
trophysical observations as limits on them may provide 
useful guidance for devel oping g rand unification theories, 
see, e.g., refs. [ill,!!!!!!^,©!!!,!^^,!!!,!!!,!!^. These 
studies have estimated various upper limits on the a. In 
the range of a = 10"^° - 10^^ and A = 10^*^ - 10"" m, 
there is a clear trend of increased strength a at shorter 
length A. What wc have constrained is the value of g^ / p^ 
or equivalcntly the \a\\'^ from the pressure necessary to 
support both static neutron stars and the fastest pulsars. 
While we expect that the range parameter A has to be 
much larger (smaller) than the radii of finite nuclei (neu- 
tron stars), we can not set separate constraints on the 
values of a and A. Compared to other efforts to con- 
strain the a and A, our study here is unique in that the 
estimated minimum value of g^ j p^ is a lower limit satis- 
fying all known constraints from both terrestrial nuclear 
experiments and observations of global properties of neu- 
tron stars. Moreover, very interestingly, our estimated 
range of g^ / p^ overlaps well with the upper limits esti- 
mated from analyzing the neutron-proton and ncutron- 
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10- 



lead scatte ring data in the range of A w 10 

Iz3,[2l|,[2l,[2i. 

In summary, neutron stars arc a natural testing ground 
of grand unification theories of fundamental forces. Con- 
sidering the possible violation of the ISL of gravity, the 
stability and observed properties of NSs can not rule out 
super-soft symmetry energies at supra-saturation den- 
sities. Given the uncertainties and model dependence 
involved in extracting information about the EOS and 
symmetry energy from heavy-ion reactions, it is very im- 
portant to test the possible super-soft symmetry energy 
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at supra-saturation densities using several obscrvablcs 
simultaneously from independent experiments analyzed 
using different models. If confirmed, it may point to- 
wards a violation of the ISL in neutron stars. 
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